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ABSTRACT 

Using high resolution cosmological A'^-body simulations, we investigate the survival 
of dark matter satellites falling into larger haloes. Satellites preserve their identity 
for some time after merging. We compute their loss of mass, energy and angular mo- 
mentum as dynamical friction, tidal forces and collisions with other satellites dissolve 
them. We also analyse the evolution of their internal structure. Satellites with less 
than a few per cent the mass of the main halo may survive for several billion years, 
whereas larger satellites rapidly sink into the center of the main halo potential well 
and lose their identity. Penetrating encounters between satellites are frequent and may 
lead to significant mass loss and disruption. Only a minor fraction of cluster mass (10 
per cent on average) is bound to substructure at most redshifts of interest. We dis- 
cuss the application of these results to the survival and extent of dark matter haloes 
associated with cluster galaxies, and to interactions between galaxies in clusters. We 
find that « 35 — 40 per cent of galaxy dark matter haloes are disrupted by the present 
time. The fraction of satellites undergoing close encounters is similar to the fraction 
of interacting or merging galaxies in clusters at moderate redshift. 
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1 INTRODUCTION 

Kinematics and dynamics of galaxies within clusters are fun- 
damental in understanding galaxy formation and evolution. 
In models of hierarchical clustering this issue is closely re- 
lated to the coUisionless dynamics of dark matter haloes 
hosting individual galaxies. The fate of these haloes, after 
merging with other haloes, depends on their relative mass, 
velocity, orbit parameters, and their internal structure in 
an interconnected way. The complexity of this problem has 
led to different approximations describing the interaction in 
terms of distinct physical processes treatable analytically. 
These processes include dynamical friction (e.g. White 1976; 
Binney & Tremaine 1987), tidal forces (Mamon 1993 and 
references therein) and resonant orbit coupling (Weinberg 
1994). 

The analytic approach provides a way to study the evo- 
lution of interacting haloes, and to estimate the survival time 
of satellites falling into the potential well of larger systems 
(Spitzer 1958; White & Rees 1978). These predictions have 
been tested with well-designed numerical models, where all 
the different parameters of the problem are under control; 
in this case, the accuracy and limitations of each approxi- 
mation can be assessed (e.g. Aguilar & White 1985; Moore, 
Katz & Lake 1996). 

However, an overall evaluation of the dynamical evo- 



lution and survival of haloes in a larger system and in a 
fully cosmological context has been investigated only re- 
cently (Klypin, Gottlober & Kravtsov 1997; Ghigna et al. in 
preparation). The evolution of clustering in a specific cosmo- 
logical model naturally provides the choice and combination 
of the different parameters: the mass distribution of haloes, 
their merging rates, orbital parameters, internal density and 
velocity profiles. 

The analysis of surviving substructure in a cosmological 
context also is of great interest since recent semi-analytical 
models of galaxy formation require a self-consistent recipe 
for the efficiency of galaxy merging (Kauffmann, White & 
Guiderdoni 1993; Baugh, Cole & Frenk 1996). 

Unfortunately, a cosmological approach to this prob- 
lem has two main inconveniences. First, the required nu- 
merical resolution, both in mass and force, is usually inad- 
equate. Secondly, the ability to follow in detail the merging 
history and fate of each satellite requires a non negligible 
amount of work and computer resources. This problem be- 
comes more serious for simulations with millions of particles, 
which would be more suited for this study. 

In the present paper, we overcome these problems by us- 
ing high resolution A'^-body simulations of individual galaxy 
clusters extracted from lower resolution cosmological simu- 
lations, as described in Tormen, Bouchet & White (1997). 



© 0000 RAS 



2 G.Tormen, A.Diaferio and D.Syer 



Our simulations have both a manageable size and sufficient 
resolution to address the survival of substructure accreted 
by galaxy clusters. The present investigation complements 
the work presented in Tormen (1997), where we studied the 
merging history of all the dark matter haloes formed in the 
simulations, and followed the haloes until they were accreted 
by the main progenitors of the final clusters. 

The plan of the paper is as follows. Section ^ briefly 
describes the simulations and the algorithm used to identify 
dark matter haloes. In Section ^ we recall the physical pro- 
cesses of interest for this problem and show a few examples of 
evolution of substructure. Section ^ deals with the evolution 
of halo orbits inside the cluster, and in Section ^ we study 
the global mass loss for the accreted satellites, and deflne 
the corresponding survival times. Section |^ studies close en- 
counters between satellites within the cluster, and estimates 
their relevance to the disruption of substructure. Section ^ 
describes the evolution of the satellite internal structure. In 
Section |^ we measure the fraction of cluster mass bound to 
substructure, and in Section ^ we discuss a few astrophysi- 
cal applications of the present results. We flnally summarize 
our results in Section 



2 METHOD 

2.1 The simulations 

We use the nine TV-body simulations described in Tormen 
et al. (1997). Each simulation models the formation of an 
individual galaxy cluster in an Einstein-de Sitter universe, 
with an initial density perturbation power spectrum P{k) oc 

. More massive particles model the external tidal field 
(see Tormen et al. 1997 for details). The Hubble constant 
is Hq = 50 km s~^ Mpc~^. Each final halo is resolved by 
~ 20000 dark matter particles within its virial radius 
each particle with mass rrip ~ 5 X 10^° M0. The effective 
force resolution is ~ O.Oliit,. An example of cluster at 2 = 
is shown in Figure |l]. 

The slope of P(k) is chosen to mimic a standard CDM 
model on the scale of galaxy clusters. However, on galactic 
scales our spectrum is shallower than the CDM spectrum, 
and this results in much more small-scale structure than one 
expects in a corresponding CDM simulation with compara- 
ble resolution. This feature creates a larger statistics for the 
substructure. As the orbital parameters of satellites accreted 
by a cluster depend very little on the cosmology (Tormen, 
Frenk & White, in preparation), our results will also hold 
for other cosmological models. 



2.2 Identification of substructure 

We identify dark matter satellites using an overdensity crite- 
rion. Each satellite groups all the particles within a sphere 
of mean density contrast p/pb = 178 (with pi, the mean 
background density), centred on the local minimum of the 
potential energy (see Tormen 1997 for details). The corre- 
sponding radius of the satellite is its virial radius 

We wish to study the fate of haloes merging with each 
other. To simplify this task, we focus our study on satellites 
accreted directly onto the main progenitor of the main halo 
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Figure 1. One of the nine clusters, at 2 = 0. Only particles 
within the virial radius are plotted. The tick spacing is 200 kpc. 
Notice the amount of resolved substructure. 



in each simulation. The main progenitor of a halo is cho- 
sen following its merging history back in time, and selecting 
the main branch at each split of the merging tree (see e.g. 
Lacey and Cole 1993). For each simulation, this procedure 
provides a catalogue of satellites which will merge with the 
main progenitor halo at different times. 

The identification time tid of a satellite is defined as the 
last time output before the satellite accretes onto the main 
halo progenitor. The merging time tmer of a satellite is the 
time when the satellite first crosses the virial radius of the 
major halo. Consecutive outputs of our simulations are sep- 
arated by an interval At ~ 1.6 Gyr. Therefore, the merging 
time for a satellite can only be bracketed between tid and 
tid + At. We assign to each satellite a nominal merging time 

tmer = tid ~t~ 0.5At. 

After a satellite is accreted by a larger halo, its orbit and 
its internal structure are perturbed by the interaction with 
the main halo and with other substructure, and the satellite 
may eventually dissolve. To keep track of its motion after 
merging, we consider two different masses associated with 
the satellite: (i) the fract ion of its initial mass which remains 
self-bound (Section 5.2), and (ii) the fraction of its initial 
mass within its tidal radius (Section ^^) . These definitions 
identify the satellite as a separate dynamical entity within 
the main halo. 

The simulations were run for a Hubble time, i.e. ~ 13 
Gyr for our cosmological model. Because we consider satel- 
lites which merged between z ~ 2.8 and z ~ 0.07, we can 
follow them within the cluster for up to ~ 11 Gyr, depend- 
ing on their merging time. We limit the present analysis to 
satellites with A'^ > 30 particles within r„ at tmer, in order 
to reduce the importance of numerical effects on the results. 
The complete set of nine simulations yields 461 satellites 
satisfying this constraint. 
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Figure 2. Snapshots of three mergers. The circle denotes the virial radius of the parent halo. Only the particles initially within the 
virial radius of the satellite are plotted in each frame. To give an idea of our resolution, the satellite in the third example is formed by 
68 particles, of which 51 are still self-bound in the final output. 



3 TIDAL FORCES AND DYNAMICAL 
FRICTION 

In this section we review qualitatively the physical processes 
dissolving the substructure. 

Dynamical friction drives satellites towards energy 
equipartition with the smooth distribution of the cluster 
particles. Satellites are more massive than single particles, 
hence they are slowed down, their orbits shrink and become 
more circular as the satellites lose energy and angular mo- 
mentum. Eventually the satellites are deposited in the center 
of the cluster potential well. 

Two different types of tidal forces must be considered 
during the infall of satellites onto a cluster: (1) global tides 
caused by the interaction with the main halo, and (2) tides 
caused by encounters between satellites. Global tides unbind 
mass from the satellite in favour of the deeper potential well 
of the main halo. For almost circular orbits this effect is 
relatively mild, and only the mass in the outskirts of the halo 
is stripped off. Conversely, satellites on very eccentric orbits 
pass through the cluster core, and are completely evaporated 
by tidal shocks. Besides global tides, close encounters with 
other satellites convert orbital energy into internal energy, 
causing collisional stripping and eventually disruption. 

Dynamical friction and tides are in action at the same 
time; however they have somewhat different effects on a 
satellite. Dynamical friction mostly influences the orbit of 
a satellite, by always reducing the orbital energy and angu- 
lar momentum, making satellites more susceptible to global 



tides. Close encounters between satellites within the clus- 
ter change both the satellite orbit and its internal structure. 
Orbits are modified randomly, with no predictable net ef- 
fect, and a satellite may either lose part of its mass or even 
capture mass from the perturber. Global tides mainly affect 
the internal structure of satellites. 

Figure ^ shows three examples of mergers which roughly 
correspond to the extreme behaviors found in the simula- 
tions. They belong to the same simulation and correspond 
to three satellites which merged with the main cluster at 
the same redshift z = 0.9. In the first row the satellite and 
the main halo have mass ratio 0.5. For massive mergers like 
this one, dynamical friction is very efficient, and the satel- 
lite is driven to the center of the main halo in a very short 
time, almost regardless of the satellite initial orbit. The ac- 
creted satellite is also heated by the tidal shock caused by 
the potential of the main halo and possibly by close en- 
counters with other satellites. However, due to its relatively 
large initial size, a significant part of its mass is confined 
within its original virial radius r^. In the second and third 
row the mergers have a much smaller mass ratio, of the or- 
der of 0.01. For this mass ratio, dynamical friction time is 
longer than a Hubble time, and tidal forces essentially deter- 
mine the disruption or survival of a satellite. In the second 
example the satellite experiences a close encounter with an- 
other satellite inside the main cluster, and evaporation is 
immediate. In the third example, the orbit avoids both the 
cluster core and collisions with other satellites, allowing a 
much longer survival time. Therefore, the disruption or sur- 
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Figure 3. Evolution of orbits: distance, angular momentum, and 
circularity. Solid, dotted, short-dashed and long-dashed lines are 
for m„/M„ e (0.00,0.01], (0.01,0.05], (0.05,0.20], (0.20,1.00], re- 
spectively, where rUv/Mv is the ratio between the satellite mass 
and the main halo mass at merging time. The light dotted line is 
at half of the initial distance and angular momentum. 



vival of a satellite depends on the balance between its initial 
orbit, close encounters with other haloes, and the efficiency 
of dynamical friction. 



4 EVOLUTION OF ORBITS 
4.1 Orbits of satellites 

In this section we study the systematic orbital decay experi- 
enced by satellites after crossing the virial radius of the main 
cluster. We consider all the time outputs after the merging 
time tmer and calculate the main orbital parameters for each 
r = t — truer - Figurc |^ shows the results. To show the depen- 
dence of the orbital decay on the satellite mass, we consider 
four subsets, according to the ratio mv{tmer) / M^(tmer) be- 
tween the satellite mass mv{tmer) and the mass Mv{tmer) of 



the main cluster at merging time. We then calculate the me- 
dian value of the parameters over the satellite distribution 
for each mass bin and for each r. 

The top panel of Figure ^ displays the evolution of the 
mean distance R{t) of a satellite from the cluster center. 
Note that distances at a given time are in units of the virial 
radius of the main halo at that time. Therefore the median 
distance decreases both because of dynamical friction and 
because of the increasing size of the main halo. Interpreting 
the figure with this caveat, we observe that the orbital de- 
cay is faster for more massive satellites, as expected from dy- 
namical friction. The time required to halve the mean orbital 
distance is roughly 5 Gyr for satellites with m„/M„ < 0.01, 
4 Gyr for 0.01 < m„/M„ < 0.05, and less than 1 Gyr for 
niv/Mv > 0.05. Due to the mass cutoff at 30 particles, no 
satellite was accreted 11 Gyr before the final time in the 
smallest mass bin. Therefore the solid curve is less extended 
than the others. It does not mean that the distance drops 
to zero at r ~ 11 Gyr. 

The second panel shows the evolution of the orbital 
angular momentum. This quantity is almost conserved for 
satellites with m^/A/^ < 0.05, while the halving time is 4 
Gyr for 0.05 < m„/A/„ < 0.20, and less than 1 Gyr for 
rUy/My > 0.2. 

The bottom panel shows the evolution of the orbital cir- 
cularity e, defined as the ratio between the orbital angular 
momentum and the angular momentum of a circular orbit 
with the same energy. This parameter is a measure of the or- 
bit shape: values of e ~ correspond to almost radial orbits, 
while circular orbits have e = 1. The energy equipartition 
caused by dynamical friction leads to orbital circularization, 
as indeed observed in the figure. However, the dependence 
on the satellite mass is weaker for circularization than for 
orbital decay and angular momentum loss. Only the most 
massive satellites reach complete circularization in a Hubble 
time. 



4.2 Comparison with theory 

Here we compare the actual orbits of the satellites in the 
Ai'-body simulation with a theoretical prediction based on 
a simple local (Chandrasekhar 1943) prescription for dy- 
namical friction. The theoretical orbits were calculated with 
time- dependent satellite masses. The orbit of each satellite 
was integrated in the time-dependent potential of its parent 
cluster as measured from the simulations. For this purpose, 
the potential was derived from the spherically averaged mass 
profile and its centre of mass was fixed. For the dynamical 
friction we applied a frictional force 



Ap{t, r)msat{t) 



(1) 



where p{t, r) is the local halo density at the position of the 
centre of mass of the satellite, Wsat is the instantaneous ve- 
locity of the satellite, and s is defined by 



s = max\vsat,Vc], 
where 

wc = {M{R)/Ry^^ 



(2) 



(3) 



is the circular velocity in the halo at the satellite position. 
Equation (^) is a crude approximation to the usual local 
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bins (the same mass bins as defined for Figure ^). In all 
cases the theoretical prediction is rather good for the two 
smallest mass bins, and fails for the large mass satellites as 
they approach the centre. This is not surprising since the as- 
sumptions behind the simple orbit model break down when 
M{R) is comparable with rrisat- The most drastic effect of 
the breakdown of the assumptions is the plunging of the the- 
oretical orbits of the large satellites towards the centre of the 
parent cluster, with J(r) and e(r) tending rapidly to zero. 
Even for these most massive satellites, however, the evolu- 
tion up to ~ 2 Gyr of R{t) and (to a lesser extent) J(r) 
is well enough reproduced by the theoretical calculations. 
Thus, the decay times can be accurately predicted. 



5 SURVIVAL TIMES 

In this section we consider the global disruption of substruc- 
ture. One way to define the survival time Tsur of a satellite 
of mass m is through the instantaneous mass loss: 



m(r) m(r) 
m(r) Am(r) 



At 



(4) 



where Am(r) — m(r+At) — m(r). A satellite followed inside 
the cluster for n time outputs will contribute n — 1 data 
points to the distribution of Tsur- We implicitly assume that 
these points are statistically independent from each other. 
Another sensible definition of survival time Taur is simply 
the time taken by a satellite to completely lose its initial 
mass m„ = m{rv;tmer)- We will use both definitions below. 

In this section we will measure the mass of a satellite in 
two different ways: (1) the mass mt(r ) = m{rt; r) within the 
satellite tidal radius rt, (Section 5.1), and (2) the satellite 
self-bound mass msbif), (Section p.2|). Unlike mt{r), the 
self-bound mass msb{T) contains irfformation on both the 
positions and velocities of the particles. For both mass mea- 
sures, we only consider the particles within the satellite's 
virial radius r„ at tmer- 



Figure 4. Evolution of orbits compared with theory assumine 
a simple local prescription for dynamical friction: as in Figure H 
distance, angular momentum, and circularity. Each panel shows 
the ratio of the theoretical prediction to the quantity measured 
in the simulation, in each case taking the median value over all 
satellites at a given time. Solid, dotted, short-dashed and long- 
dashed lines are for m„ /M„ e (0.00,0.01], (0.01,0.05], (0.05,0.20], 

(0.20, 1.00], respectively. 



Chandrasekhar formula for dynamical friction in an isotropic 
system (e.g. Binney & Tremaine 1987). Near the centre of 
the parent halo Vc is softened by replacing R with \/ R^ + e^, 
where e ~ 2kpc. The Coulomb logarithm A was set to a con- 
stant (= 8). We tried two prescriptions for the satellite mass 
rUsat- 1) the 'virial' mass mv{t), more exactly the mass of 
satellite particles within the original virial radius at the iden- 
tification time; and 2) the self-bound mass m(,(t), defined in 
Section 5.2. The former was found to give the best results (a 
similar conclusion was reached by Navarro, Frenk & White 
1995), and was used for the results described. 

The results are shown in Figure ]A Here we show the 



ratio of the theoretically predicted to the measured quanti- 
ties in terms of the median over all satellites in each mass 



5.1 Tidal radius 

Consider a satellite in a circular orbit at distance R from the 
cluster center. The orbiting satellite is tidally truncated at 
some radius rt, loosely speaking where the differential tidal 
force of the cluster is equal to the gravitational attraction 
of the satellite: 



dFud{R) ■ 



dF 
'dR 



dR = 



Gm 



(5) 



Assuming that the satellite mass m and its radius rt are 
negligible compared with the cluster mass and with their 
relative distance R, i.e. m <^ M{R), and rt = dR <^ R, we 
readily obtain 



rt 



R 



(2 - d\nM/d\nR)M{R) 



1/3 



(6) 



Therefore the tidal radius is such that the mean density of 
the satellite within rt is of the order of the mean density of 
the main halo within R. This definition captures the essence 
of the natural definition of rt, defined as the distance of the 
center of mass of the satellite from the saddle point of the 
potential of the total system. 
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Figure 5. Left panels: evolution of tidal radii. Bold solid line is the median evolution; shaded areas include first and third quartiles. 
Right panels: differential (histograms) and cumulative (solid lines) distributions of the survival times defined in Equation (^) for the 
same mass bins of the left panels. We remind that is the virial mass of satellites at merging time tmer, and My is the virial mass of 
the main cluster progenitor always at tmer- 



For non circular orbits, the most common situation in 
real life, the satellite tidal radius is rather ill-defined (e.g. 
King 1962; Binney and Tremaine 1987; Mamon 1987; Ma- 
mon 1993). In this case, one usually keeps Equation (^, tak- 
ing for R the pericentric distance of the satellite orbit. We 
adopt this procedure. At each time r, rt is the solution of the 
equation psat{rt\ r) = pmain{rp; r - At), where psat{rt; r) is 
the satellite mean density at rt, and Pmain{rp;T — At) is 
the mean density of the main halo at the pericenter Vp of 
the orbit the satellite had at the previous time output. We 
compute density profiles for the satellites at each time con- 
sidering only the particles within at tmer- 

We then define the survived mass of a satellite as the 
mass m{rt; r) within its tidal radius. With this mass. Equa- 
tion (^) provides our first definition of survival time Tsur of 
a satellite. Note that the tidal radius is usually associated 
with the effect of global tides on a static satellite profile. We 
use the actual density profile of the satellite, which changes 
in time accordingly to all kind of interactions experienced by 
the satellite. Therefore, the corresponding survival times are 
a measure of the times taken by a satellite to be destroyed 



by all processes together: i.e. cluster tides, close encounters 
and dynamical friction. 

Figure ^ shows the evolution of tidal radii rt and the 
distribution of survival times for the accreted satellites. We 
bin satellites into four mass bins as in Figure ^ As one could 
expect, tidal radii of small satellites vanish more slowly than 
tidal radii of massive satellites. Notice that in all mass bins 
rt drops to 30 - 50 per cent of its initial value as soon as the 
satellite enters the main cluster and finds itself embedded 
in the denser environment. Because of the sudden change in 
the surrounding density, particles with large kinetic energy 
escape the satellite. The dense core survives, and afterwards, 
rt declines to zero more gently, with an intermediate con- 
stant phase for the less massive satellites. 

The corresponding survival times are displayed in the 
right panels. The distributions are plotted up to a Hubble 
time, but they have a long tail further on the right, as indi- 
cated by the cumulative distributions. Median survival times 
for the four mass bins are, from low to high mass, Tsur = 7 
Gyr, 5.5 Gyr, 2.4 Gyr and less than 1.5 Gyr. The fraction 
of survival times above 13 Gyr is 36 per cent, 30 per cent, 
10 per cent and 3 per cent from low to high mass. 
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5.2 Self-bound mass 

Another measure of the mass associated with substructure is 
the fraction msb{T)/mv{0) of the initial satelhte mass which 
remains gravitationally self-bound. The self-bound mass of 
a satellite was defined as follows. 

(i) At any given time we consider only the particles which 
composed the satellite at its identification time. 

(ii) We estimate the total energy of each particle summing 
the potential energy due to the distribution of these parti- 
cles and the kinetic energy calculated in a reference frame 
moving with the average velocity of all particles within the 
original virial radius rv of the satellite; the reference frame 
is centered on the position found by the moving center tech- 
nique (see Tormen et al. 1997 for details). 

(iii) We remove all particles with positive total energy, 
and calculate a new center of mass and average velocity for 
the distribution of the remaining particles. 

(iv) we calculate new total energies using the new set of 
particles and the new center of mass and velocity. 

We iterate the last two steps until the number of particles 
with negative energy is constant. This final value gives the 
self-bound mass of the satellite ms6(r). 

In Figure ^ we compare the evolution of the self-bound 
mass with the evolution of the mass within the tidal radius. 
Bold lines show median values, while shaded areas show the 
first and third quartile range of the distribution at each time 
T. The trend of the two estimates is similar, but the mass 
within rt is always smaller than the self-bound mass. 

It is not surprising that the two definitions do not agree. 
In fact, they use rather different criteria to identify the satel- 
lites. The mass within the tidal radius is defined by the spa- 
tial distribution of particles alone, Eq. the self-bound 
mass also uses information coming from the particle veloci- 
ties. Because particles in the outskirt of a satellite will take 
roughly a free-fall time to be accelerated to the cluster veloc- 
ity dispersion, they will stick to the satellite for some time 
after they fall outside the tidal radius, and will be counted 
in the self-bound mass. Thus, the estimate based on tidal 
radii is actually a conservatively small measure of the mass 
associated with the satellite. 

We use the evolution of the mass within rt as our sec- 
ond estimate of survival times. Figure ^ shows that median 
survival times, defined by the condition: m(rt; Tsur) = 0, are 
Tsur = 7.5, 7.5, 4 and 2.5 Gyr from small to large masses. 
These survival times are slightly larger than the survival 
times derived with Equation (Q). However, they are consis- 
tent with those survival times, because both distributions 
have large scatter. 

The longer survival time of small satellites is due to the 
combination of two effects: (1) smaller satellites are more 
compact (Tormen 1997; see also Section ^; (2) dynamical 
friction is less effective on smaller satellites; in fact, smaller 
satellites have larger distances from the cluster center (Fig- 
ure and suffer a weaker global tide. 

These results suggest that the high force and mass res- 
olution of our simulations overcome, at least partially, the 
overmerging problem which is common to dissipationless A'^- 
body simulations (Klypin et al. 1997). Moreover, survival of 
galaxy size haloes do not necessarily need dissipative simu- 
lations as usually believed (e.g. Summers et al. 1995). For 



example, satellites with mass ratio below 0.01, roughly corre- 
sponding to galaxy size haloes falling into a forming cluster, 
have median survival time 7-7.5 Gyr. Thus, half of the satel- 
lites of this size, accreted by a cluster at 2: ~ 0.6 — 0.8, are 
still "safely" orbiting within the cluster potential at 2: = 0. 
One example is the third satellite in Figure ^ Defining sur- 
vival times based on self-bound masses would result in even 
longer survival. 

Finally, we find that survival times do not depend on 
the initial orbit of satellites. Radial orbits drive satellites 
through the cluster core where tidal forces dissolve them. 
However, radial orbits are rare (Tormen 1997), and our 
statistics is too poor to investigate this issue satisfactorily. 



6 ENCOUNTERS BETWEEN SATELLITES 

In this section we investigate the importance of close en- 
counters between satellites orbiting within the parent clus- 
ter. High speed close encounters have been advocated as a 
major mechanism for the morphological evolution of galaxies 
in clusters (Moore et al. 1996). The mass loss in satellite- 
satellite encounters depends on different parameters. Here 
we restrict our analysis to two quantities: the relative dis- 
tance b of the encounter, and the ratio 7 between the relative 
tangential velocity dvtg of the two satellites and the inter- 
nal velocity dispersion of the perturber. Closer (smaller 
6) and slower (smaller 7) encounters are more effective in 
disrupting the colliding satellites. 

Operationally, we look for encounters between satellites 
in the following way. Consider a satellite i accreted by the 
main cluster at time ti . The other satellites in the main halo 
are the 'perturbers' of satellite i. At each time output t > ti, 
we compute the distance of satellite i from all perturbers. 
The dimensionless relative distance between satellite i and 
perturber j at time t is b{i,j;t) = d{i, j;t) /[r{i) + r{j)\, 
where d{i,j;t) is the distance between the two satellites at 
that time, and r{i),r{j) are the virial radii of the satellites 
at their respective merging times. The dimensionless rela- 
tive velocity of the encounter ^{i,j;t) = dvtg{i,j;t)/av{j) 
is similarly defined. The fraction of self-bound mass lost by 
the satellite between the time of the encounter and the next 
time output is 

Am{i) _ msb{i;t) — m.sbii;t + At) 
m(i) msb{i;t) 

We then consider the distribution of these quantities for all 
satellites, at all times t > tmer- With the present definitions, 
a satellite may provide more than one data point to the 
distributions. We make the assumption that each satellites 
provides statistically independent data points in different 
time outputs, as was assumed in Section ^. We restrict our 
search to: 

(i) pairs of satellites retaining a non negligible fraction of 
their original mass at the time of the encounter: if fsb{t) = 
nT-sbit) /niv{tmer) IS the fraction of the initial satellite mass 
which remains self-bound at time t, we require that both 
fsbii;t) and fsbU;t) > 0.2. 

(ii) satellites found within the virial radius of the main 
halo: R{T)/R^{t) < 1. 

The first requirement excludes fake encounters between 
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Figure 6. Survived mass fractions after merging. Solid and dashed lines are the median evolutions of the self-bound mass and the mass 
within the tidal radius, respectively. Shaded areas include first and third quartiles. Mass bins are as in Figure ^ 



haloes which have aheady largely dissolved. The second re- 
quirement excludes time outputs when a satellites has left 
the main halo after a first identification. We also exclude 
satellites identified at the penultimate time output of the 
simulation, because such satellites are found within the clus- 
ter in just one output (the last one), and so we cannot cal- 
culate their mass loss. 

The top panel of Figure |^ shows the differential and 
cumulative distributions of the number of encounters per 
satellite. We asked each satellite: "how many encounters did 
you experience which satisfy the requirements listed above?" 
The answer is illustrated by the dotted curves, which refer 
to any value of b. We also asked: "how many encounters with 
6 < 1 (referred to as penetrating encounters) did you experi- 
ence which satisfy the same requirements?". The answer is 
given by the solid curves. The dotted cumulative distribu- 
tion shows that 25 per cent of the satellites accreted by the 
cluster have no encounter at all that satisfy the conditions 
listed above, while the remaining 75 per cent have 1 to 4 en- 
counters. More interesting, the solid cumulative distribution 
shows that almost 60 per cent of all satellites have at least 
one penetrating encounter. We repeated the same census 
for satellites with m^/My < 0.01, representing galaxy-size 
haloes merging with a forming cluster. Almost 85 per cent 



of these haloes have at least one encounter, and 55 per cent 
have at least one encounter with b < 1. In the other two 
panels we show the differential and cumulative distributions 
of relative distance and relative velocity for satellites hav- 
ing at least one encounter. The central panel illustrates that 
penetrating encounters are very common within the cluster: 
over 60 per cent of the encounters have fe < 1. The bottom 
panels shows that the relative velocity of the encounter is 
generally much larger than the internal velocity dispersion 
of the satellites. In fact, the satellites orbit inside the par- 
ent cluster at a speed of the order of the cluster velocity 
dispersion, which is much larger than the internal velocity 
dispersions of the satellites themselves. As a consequence, 
only a few per cent of the encounters have 7 of order unity. 
This high relative speed justifies the approximation of ignor- 
ing mass capture during the encounters, which we implicitly 
make when we assigning particles to a satellite. 

To disentangle the mass loss due to low speed, close 
encounters, from the mass loss due to the tidal forces of the 
main halo, in Figure |^ we consider mass loss versus b (top 
panel) and mass loss versus relative velocity (bottom panel), 
for encounters occurring at different distances R{t) / Rv from 
the center of the main halo. The top panel shows that, for 
any fixed b <; 0.2, the mass loss is roughly independent of 
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Figure 8. Top panel: median mass loss versus relative distance 
of the encounters b. Different symbols refer to encounters oc- 
curring at different distances from the center of the main halo: 
< R{t)/Rv < 0.4 (solid squares), 0.4 < R(t)/Rv < 0.7 (open 
squares) and 0.7 < R{t)/ R^j < 1 (open circles). In each curve the 
data are binned so that each point represents the same fraction of 
the total encounters. Error bars are ilcr of the median. Bottom 
panel: median mass loss versus relative velocity of the encounter. 
Symbols and data binning as in the top panel. 



Figure 7. Parameters in satellite-satellite encounters. In each 
panel the thin and thick histograms indicate the differential and 
cumulative distribution respectively. Top panel: total number of 
encounters (dotted) and of penetrating encounters (solid) expe- 
rienced by satellites within the main halo. Middle panel: Di- 
mensionless relative distance of the encounter. Bottom panel: 
Relative tangential velocity of the encounter, in units of the in- 
ternal velocity dispersion of the perturber. 



this distance: encounters in the central region of the main 
halo (0 < R{t)/Rv < 0.4: solid squares) and those taking 
place in outer regions (0.4 < _R(r)/_R„ < 0.7: open squares, 
and 0.7 < R{t)/R^ < 1: open circles) cause the same mass 
loss, within the statistical noise. 

Notice the different range of b covered by the curves: the 
solid symbols are concentrated at fe < 1, while encounters 
happening at larger distances extend further to the right. 
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In fact, the number density of satellites grows towards the 
center of the main halo, increasing the chance of close en- 
counters. The solid square at Am/m ~ 0.9 corresponds to 
head-on encounters in the inner part of the main halo, which 
lead to total mass disruption. 

For 6^1 the curves are consistent with a roughly con- 
stant mass loss: Am/m ~ 0.15. If we interpret this as the 
average mass loss due to global tides only, then the extra 
mass loss is due to the encounter. The fact that mass loss 
increases only for penetrating encounters {b < 1) is consis- 
tent with this explanation. 

The lower panel shows that the mass loss is almost in- 
dependent of the value of the relative velocity 7, and of the 
distance R{t)/Rv from the center of the main halo. The ex- 
ception is again the leftmost solid square, which corresponds 
to slow encounters in the core of the main halo. Encounters 
in the core of the main halo are both the closest and the 
slowest, hence the most disruptive. In fact, (1) the velocity 
dispersion and circular velocity of the main haloes decrease 
towards the center for r/i?„ ^ 0.3 (as shown by the radial 
profiles in Tormen et al. 1997), and (2) the number density 
of satellites increases towards the center. 

On the whole, Figure |^ indicates that the mass loss dur- 
ing encounters between satellites depends mainly on the rel- 
ative distance of the encounter, b, and very little on the rel- 
ative velocity, because slow encounters are rare, while close 
encounters are frequent. Therefore, in the rest of this section 
we will concentrate our study on the b parameter. 

Since we found that a constant mass loss Am/m = 0.15 
is a reasonable zeroth order description of the effect of global 
tides, we can estimate the mass loss due only to collisions 
by subtracting this value from that originally measured: 

f^V = ^-0.15. (8) 

\ m J m 

The distribution of this quantity is equivalent to the disrup- 
tion time associated with the encounters: 

'mir) ( m(r) Y 

Taou = --rj-T = - -r — i—r At. (9) 
m(r) \ Am(T) J 

The top panel of Figure |^ shows the differential and 
cumulative distributions of TcoU for penetrating encounters 
(b < 1: solid lines) and for encounters with 6 > 1 (dotted 
lines). For this distribution we used all the data of Figure ^ 
excluding those contributing to the leftmost solid symbol. 
The median disruption time for 6 < 1 is Tcou — 11 Gyr, 
corresponding to a median mass loss of (Am/m)' ~ 0.15. 
The percentage of satellites dissolved by collisions over a 
Hubble time is 52 per cent. On the other hand, encounters 
with 6 > 1 have negligible mass loss, and median Tcoh 3> 
t Hubble, a result consistent with the data shown in the top 
panel of Figure ^ 

The bottom panel of Figure ^ shows the same statis- 
tic for satellites with m„/M„ < 0.01, which correspond to 
galaxy-size haloes falling onto a cluster. In this case mass 
losses due to penetrating encounters are smaller, with me- 
dian value Am/m = 0.05, corresponding to a median dis- 
ruption time t ~ 30 Gyr, with 42 per cent of satellites being 
dissolved over a Hubble time. 

We can use the data in the first panel of Figure ^ on the 
frequency of penetrating encounters in satellites, to quantify 
the statistical significance of encounters. Since 60 per cent 
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Figure 9. Differential (thin histograms) and cumulative (thick 
curves) distributions for the disruption time t^oH from Eq. (^ 
due to close encounters between satellites. Top panel: Distri- 
bution for satellites of any mass. Solid and dotted lines are for 
b £ [0,1), [1,00) respectively. Bottom panel: Same as in top 
panel, for satellites in the mass range niy/M^ < 0.01. Values of 
P12 in each panel are the Kolmogorov-Smirnov significance levels 
for the difference between the solid and dotted distribution. 



of all satellites have one or more penetrating encounters, 
from the top panel of Figure ^ we can say that at least 
30 per cent of all satellites are dissolved in 11 Gyr or less 
by penetrating encounters alone, and a comparable fraction 
is dissolved over a Hubble time. Similarly, 23 per cent of 
satellites with m^/Mv < 0.01 are dissolved by penetrating 
encounters in less than a Hubble time. These numbers are 
an estimate of mass loss and disruption time due only to 
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penetrating encounters. Total mass losses and total survival 
times are those presented in Section ^. 

Finally, we found that these distributions are fairly ro- 
bust to variations of parameters. We obtain essentially the 
same results if we change the requirements on the fraction 
of initial satellite mass which must be self-bound at the time 
of the encounter, or if we only consider encounters at larger 
distances R{t) / Rv from the center of the main halo. 



7 INTERNAL STRUCTURE: MASS AND 
VELOCITY PROFILES 

Here, we investigate the evolution of the internal structure 
of satellites. All the particles composing the satellite just 
before the merging time define the density and the velocity 
profile at each time. We can thus investigate how the internal 
structure of a satellite changes after the merging time, and 
how its mass is redistributed within the main halo. 

The internal structure evolution of satellites is different 
for small and large mass satellites. The density of the en- 
vironment surrounding a satellite sufi'ers an abrupt change 
when the satellite enters the main halo. The tidal radius 
drastically decreases (Figure ^ and particles with enough 
kinetic energy escape the satellite. This process leaves low 
mass satellites with particles having a smaller rms veloc- 
ity (Figure [ll] ), producing an overall effect of cooling and 
a slight mass concentration (Figure |l^). Escaped particles 
quickly thermalize to the main halo velocity dispersion. 

From Figure |lo| we see that small satellites are more 
concentrated than massive ones. In fact, half-mass radii 
for the initial profiles are at r/r„ = 0.2, 0.28, 0.35, 0.41 
from small to large masses. For this reason, when we con- 
sider more massive satellites, the equality between the mean 
density of satellites and that of the main halo occurs at 
lower values of r/r„. Thus, massive satellites are heated well 
within r„, and the cooling effect disappears. 

Moreover, the initial rms velocity of massive satellite 
particles is closer to the rms velocity of the main halo; thus, 
their kinetic energy gain is less dramatic. Mass profiles also 
show that the "disruption" of massive satellites is actually 
a mere infiating (see also Figure ^). In fact, the mass pro- 
files are simply shifted to larger values of r, but their shape 
remains similar. On the other hand, the external shells of 
small satellites are violently stripped off and particles are 
redistributed within the main halo. 



8 THE FRACTION OF HALO MASS IN 
SUBSTRUCTURE 

From our study it appears that only relatively small satel- 
lites can resist dynamical friction and survive in a cluster 
for a significant time. However, most of the mass forming 
a cluster comes from a few, massive objects (e.g. Tormen 
1997). Therefore, it is natural to ask what is the fraction of 
cluster mass which, at any given time, is bound to substruc- 
ture. The mass m(rt; r) within the tidal radius is a measure 
of such fraction. By definition, the mean density within rt 
is larger than the mean density of the main cluster within 
the radius corresponding to the pericentre of the satellite. 
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Figure 12. Fraction of cluster mass in substructure. Top panel: 
each curve represents the median (over the cluster sample) num- 
ber of satellites found within the main halo, with mass larger 
than the threshold indicated in the figure, as a function of red- 
shift. Bottom panel: Median values (over the cluster sample) of 
the fraction of cluster mass associated to satellites, as a function 
of redshift. Open circles, open squares and solid squares refer to 
the fraction within a sphere of radius 0.2R^, 0.5Rv and i?„ re- 
spectively. Masses are defined using tidal radii, as in Section 3.1. 



Therefore, the mass within the tidal radius always corre- 
sponds to an enhancement in the local mean density. We 
refer to this as a 'visible' structure. By this definition any 
satellite with rt > is visible. 

The top panel of Figure ^ shows, as a function of 
redshift, the median number of satellites within the main 
halo, and with mass m(rt; t)/M-u{t) above a given threshold; 
Mvit) is the virial mass of the main halo at that time. It 
shows that, on average, there are 2 satellites with tidally 
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Figure 10. Evolution of satellite mass profiles. Lines refer to different times. Time sequence is as follows: solid, dotted, short-dashed, 

dot-dashed, dot-dot-dot-dashcd, long-dashcd, bold long-dashed, solid. Note how the mass profile of the small satellites (top panels) is 
severely disturbed at later times (long-dashed, bold long-dashed, solid curves). Massive satellites (bottom right panel) merely expand 
self-similajrly. 



truncated mass larger than O.OIM^, one with mass larger 
than 0.02M„ , while there is no satellite larger than 0.05M„ . 

Note that these mass thresholds roughly correspond to satel- 
lites with radius 20, 27 and 37 per cent that of the main halo. 
Such objects should be fairly easily identified by e.g. weak 
lensing observation of the dark matter distribution in galaxy 
clusters (see e.g. Geiger & Schneider 1997). 

The mass attached to satellites is shown in the bottom 
panel. The data show that, for most redshifts of interest, 
substructure make up on average 10 per cent of the cluster 
mass within Ji„ (solid squares). The maximum value mea- 
sured is 30 per cent. In the cluster inner regions the fraction 
of mass in substructure lowers dramatically (open squares 
and open circles refer to spheres of radius 0.5-R„ and 0.2_R„, 
respectively), due to the stronger tidal forces and to the nu- 
merical resolution. This estimate is really a lower limit, as 
the mass within rt is a conservatively low estimate of the 
mass of a satellite. Nevertheless most of the cluster mass 
is not in substructure, but is smoothly distributed inside 
the cluster. Therefore the issue of survival of substructure 
in massive haloes mainly applies to small {my/ My <, 0.05) 



and compact satellites rather than to large {my/My ^ 0.05) 
ones. 



9 DISCUSSION 

9.1 The extension of dark matter haloes in cluster 
galELxies 

Studies of the kinematics of satellites of field spirals indi- 
cate that these galaxies are embedded in dark halo at least 

10 times larger than the optical radius of the galaxy (Zarit- 
sky & White 1994; Zaritsky et al. 1997). In cluster galaxies, 
such haloes are thought to be tidally-truncated by the clus- 
ter potential. Unfortunately, observations of this effect are 
very difficult. As for spiral galaxies, they are depleted of their 
gas by ram pressure in the intracluster medium soon after 
they merge, so that rotation curves cannot be measured at 
distances significantly larger than the optical radius (Gunn 
& Gott 1972, White et al. 1991). The dark haloes of ellipti- 
cals may in principle be probed by the kinematics of their 
globular clusters, but these objects are quite faint, so that 
HST observations are needed to study globulars, even in 
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Figure 11. Same as Figure |l^ for the rms velocity profiles. In the top left panel, note how particles first stream out of the satellite 
(dotted curve) and they increasingly thermalize with the main halo rms velocity later on (short-dashed, dot-dashed, dot-dot-dot-dashed, 
long-dashed, bold long-dashed curves). The rms velocity variation is abrupt, indicating that satellites are physical entities distinguishable 
from the main halo. The thermalization process is smoother with increasing satellite mass. The central part of the most massive satellites 
(bottom right panel) suffer very little cooling. Finally, note that, on average, small satellites are completely disrupted at later times (bold 
long-dashed, solid curves). 



nearby cluster galaxies. If dark haloes turn out to be small 
enough, tidal truncation could be directly observable using 
weak gravitational lensing (Geiger & Schneider 1997) or ob- 
serving a suitable tracer population to large radii. 

Lacking obser vati onal evidence, we can at least use the 



results of Section 5.1 to constrain the expected radius of 
dark matter haloes in cluster galaxies. We can associate a 
galaxy radius to each dark matter satellite using the dataset 
of Burstein et al. (1997). These authors showed that the 
dynamical properties of most stellar systems, ranging from 
globular clusters to galaxies to galaxy groups and clusters, 
obey relations similar to the Fundamental Plane of ellipti- 
cals (Dressier et al. 1987; Djorgovski & Davis 1987). The 
ensemble of these planes, termed the cosmic metaplane, is 
an expression of the virial theorem, tilted by variations in 
mass-to-light ratio. In particular, Burstein et al. present data 
for ~ 900 elliptical and spiral galaxies, for which they derive 
effective radii re (i.e. radii enclosing half of the total light 
from the galaxy) and central velocity dispersions tJc. 

Our procedure to assign re to each satellite is as follows: 



(i)we consider all satellites infalling onto our simulated clus- 
ters, and take the one dimensional rms velocity within the 
virial radius, measured at merging time, as a rough esti- 
mate of (Tc- (ii) From the sample of Burstein et al. (1997) 
we select all objects having ac within 10 per cent of the 
value found for the satellite, (in) We randomly draw one 
of these objects and assign the corresponding value of re 
to the satellite. This prescription can be applied also to a 
restricted subsample or class of observations, e.g. only the 
spirals or ellipticals in the Burstein et al. database. In such a 
case, satellites with a value of ctc not matched by any of the 
selected observations are excluded, leaving us with satellites 
in a range of tJc appropriate to the choice made. 

The result of this exercise is shown in Fig [13. In the 
top panel we plot the differential (thin histograrn) and cu- 
mulative (thick curve) distribution for the virial radius of 
satellites, in units of re, with re chosen from either elliptical 
or spiral galaxies. This distribution gives an idea of the ra- 
dius of galaxies hosted in dark matter haloes. It shows that 
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Figure 13. Radius of darlt matter satellites in units of tlie effec- 
tive radius Ve of tlic fiost-galaxy. Top panel: virial radius, mea- 
sured at merging time. Bottom panel: tidal radius, measured at 
the present time (z = 0). Details about the choice of are given 
in the main text. Thin histograms refer to differential distribu- 
tions, thick curves to cumulative distributions. Curves include all 
the nine cluster simulations. 



the median virial radius for merging satellites is rhoughly 20 
times larger than Ve. 

In the bottom panel we plot similar distributions for 
the final tidal radius of satellites (i.e. measured at 2 = 0). 
The solid distribution refers to all satellites, with re chosen 
from galaxies, galaxy groups and galaxy clusters, in order to 
cover the all range of satellite masses. The dotted distribu- 
tion refer to satellites associated to Ve chosen from elliptical 
galaxies, while the dashed distribution has chosen from 
spiral galaxies. The panel shows that ~ 35 — 40 per cent of 
galaxy-sized dark matter satellites ever accreted by a clus- 
ter have been dissolved by the present time, i.e. they have 



rtid = (dotted or dashed curves). The median tidal ra- 
dius for galaxy-sized satellites is rud — lOr-g. If we consider 
all satellites (i.e. also larger than galaxy-size), half of the 
satellites are dissolved hy z — (solid curve). 

Strictly speaking, these figures should be regarded as 
upper limits, because an increase in numerical resolution or 
the presence of a dense core of baryonic matter would prob- 
ably facilitate the survival of satellites. The relevance of the 
latter effect is however not clear, as dark matter is the dom- 
inant source of gravity already at the optical radius of the 
galaxy (e.g. Navarro, Frenk & White 1996). Taking median 
values as robust estimators of the distributions, this result 
on the whole suggests that 50 per cent of galaxies in present- 
day clusters should have a dark matter halo truncated at or 
below ~ lOre, a radius close to those radii testable by e.g. 
high precision surface photometry. Of course, some of these 
galaxies will have lost a substantial part of their dark matter 
halo and may themselves be disrupted. 



9.2 Interacting galaxies in clusters 

Observations indicate that the fraction of morphologically 
disturbed galaxies or interacting galaxies in clusters typ- 
ically increases with redshift (e.g. Lavery & Henry 1988; 
Oemler, Dressier & Butcher 1997). Due to the high rela- 
tive speed of most encounters, illustrated in Figure H, merg- 
ers between galaxies are expected to be rare. On the other 
hand, the numerical simulations of Moore et al. (1996) have 
shown that repeated high speed encounters between spirals, 
at a relative distance less than a few optical radii (e.g. ^ 50 
kpc for an L, galaxy) can seriously modify the galaxy mor- 
phology. This distance is also the typical separation of the 
interacting galaxies found by Lavery & Henry (1988). 

We can us e th e results of the analysis performed in 
Sections ^ and 9T to estimate the fraction of interacting 
galaxies in our simulations. Figure ^ shows the fraction of 
galaxy-size satellites undergoing close encounters with other 
satellites, as a function of redshift. We show results for di- 
mensionless relative distances b < 0.5, b < 0.2 and b < 0.1, 
which correspond to encounters at distances < lOre, < 4re 
and < 2re if one takes the median value rv/ve ~ 20 found in 
the previous section as representative. Encounters at such 
distances already fall in the class discussed by Moore et al. 
We see for example that ~ 20 per cent of these satellites 
are engaged in encounters with b < 0.5, while ~ 5 per cent 
are having encounters with b < 0.2, and 1 — 2 per cent en- 
counters with b < 0.1. These figures are consistent with the 
10 — 20 per cent of merging and interacting galaxies observed 
by Oemler et al. (1997) in four clusters at 2 ~ 0.4. The self- 
similarity of our simulations should ensure that such effects 
do not exhibit any trend with redshift, as it is indeed the 
case at least for z < 1. 

Our results support the claim that close encounters be- 
tween galaxies in clusters are fairly frequent. The rate of 
encounters produced by our simulations is lower than the 
one invoked by Moore et al. for galaxy harrassment, as most 
satellites in our simulations have at most one encounters at 
6 < 1 (as shown in the top panel of Figure |^ , and not sev- 
eral as required by galaxy harrassment. However, the rate 
we found should be considered as a lower limit, as in real 
clusters collisions mostly take place in the central, densest 
region of the cluster, where we do not have enough resolu- 



© 0000 RAS, MNRAS 000, 000-000 



Survival of substructure 15 



73 
CD 



CD 
-i-> 



O 
•I— I 
-i-> 

o 



0.3 



0.2 



0.1 







"1 — I — I — I — I 1 1 1 r 



b<0.5 




(l+z) 



Figure 14. Interacting satellites. The panel shows the fraction of 
satellites with m„/M„ < 0.01 experiencing encounters with other 
satellites, as a function of redshift and binned according to the 
dimensionless relative distance b. Circles refer to encounters with 
b < 0.5, open and solid squares to encounters with b < 0.2 and 
6 < 0.1 respectively. 



tion to preserve the satellite identity. Further investigation 
with higher resolution simulations is required to estimate 
this effect. 



9.3 Comparison with other work 

There exist two different approaches to studying the prop- 
erties of substructure within dark matter halos: (i) we can 
identify satellites before they merge with the main halo and 
keep track of their particle trajectories after merging; (ii) 
we can identify overdensity peaks within the main halo at 
each output time. The latter approach has the advantage 
of identifying physically motivated substructure. Moreover, 
this approach is equivalent to the observational procedure 
when it is applied to two-dimensional projections. However, 
this method is not straightforward if we are interested in the 
evolution of the substructure properties. In fact, in order to 
check the identity of density peaks identified at different 
times, we need to check whether a substantial fraction of 
the particles of a density peak at a later time belonged to a 
density peak at an earlier time (Klypin et al. 1997). 

In the present work, we have eight output times for 
each of the nine simulations. Therefore, we use the former 
approach, which is more straightforward in keeping track of 
each satellite. The relatively large number of output times is 
the most important feature of our analysis. In fact, we can 
consistently compare the satellite orbit evolution with the 
dynamical friction predictions. Moreover, we can accurately 
estimate the evolution of the satellite internal properties as 
tidal effects and collisions disrupt them. This information 
provides a quantitative definition of the survival time. 

Our analysis complements Ghigna et al. analysis. 



Ghigna et al. identify overdensities in a single halo simula- 
tion with higher space and force resolution than ours. Thus, 
they can accurately estimate the satellite density profiles, 
and the spatial distribution of the satellites within the main 
halo. On the other hand, our work has the advantage of a 
better statistics provided by the nine independent clusters 
and the use of a larger number of output times. 



10 CONCLUSIONS 

The results of this paper may be summarized as follows: 

(i) The orbital decay of haloes within galaxy clusters is 
consistent with the expectations of dynamical friction. Sub- 
structure are driven to the center of the main halo in less 
than a Hubble time if their initial mass is larger than one 
per cent of the mass of the main cluster. 

(ii) Substructure retain their identity for a significant 
fraction of the Hubble time if their initial mass is smaller 
than 5 per cent of the main cluster mass. Median sur- 
vival times, based on the mass within the tidal radius, 
are in the range [7,7.5] Gyr, [5.5,7.5] Gyr, [2.4,4] Gyr and 
[1,2.5] Gyr for mass ratios m„/M„ G (0.00,0.01], (0.01,0.05], 
(0.05,0.20], (0.20,1.00], respectively. However, the mass 
within the tidal radius is conservative, since survival times 
are more than 50 per cent longer when we consider the satel- 
lite self-bound mass. Smaller satellites have longer survival 
times for a combination of two reasons: a) they are more 
compact, and b) they are less influenced by dynamical fric- 
tion, and avoid the cluster core. 

(iii) Encounters between satellites within the cluster are 
frequent and lead to mass loss comparable to that caused 
by global tides. The mass loss is correlated with the rela- 
tive distance and almost uncorrelated with the satellite rel- 
ative velocity. In fact, slow encounters are rare, but close 
encounters are frequent. Almost 60 per cent of the satellites 
experience at least one penetrating encounter with another 
satellite before losing 80 per cent of their initially self-bound 
mass. The median mass loss per penetrating encounter is 
Am/m ~ 0.15, corresponding to a median disruption time 
Tcoii — 11 Gyr, due only to penetrating encounters. The 
same figures become Am/m ~ 0.05 and Tcou — 30 Gyr for 
galaxy-sized satellites. 

(iv) The evolution of the satellite internal structure de- 
pends on the satellite mass: smaller satellites easily loose 
their less bound particles and cool in their inner region, while 
larger satellites experience a global heating. 

(v) The fraction of cluster mass in tidally-defined sub- 
structure is 10 per cent on average within the virial radius, 
and lower in the inner parts. It therefore constitutes only a 
minor fraction of the total cluster mass. 

The application of our results to galaxies in clusters 
requires us to specify a procedure for populating the smaller 
dark matter satellites with stellar material. We discuss an 
empirical method of assigning galaxies to haloes in Section 
^ The results can be summarised as follows: 

Roughly 50 per cent of galaxies in present-day clus- 
ters should have a dark matter halo truncated at or below 
~ lOre. Some of these galaxies may themselves have been 
disrupted. 
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The fraction of satellites undergoing very close encoun- 
ters, b < 0.2 — 0.5, is similar to the fraction of interacting 
or merging galaxies in clusters at moderate redshift. Re- 
peated close encounters, as required by galaxy harrassment, 
are however very rare with the present numerical resolution. 



Zaritsky D., White S.D.M., 1994, ApJ, 435, 599 
Zaritsky D. Smith R.. Frcnk C .S.. White S.D.M., 1996, Apj sub- 
mitted, a,stro-ph/9611207 



ACKNOWLEDGEMENTS 

We thank Joerg Colberg, Riccardo Giovanelli, Martha 
Haynes, Julio Navarro and Simon White for useful com- 
ments and suggestions. Financial support for G.T. was pro- 
vided by an MPA guest postdoctoral fellowship and by 
the Training and Mobility of Researchers European Net- 
work "Galaxy Formation and Evolution". A.D. holds the 
grant ERBFMBICT-960695 of the Training and Mobility 
of Researchers program financed by the European Com- 
munity. The simulations were performed at the Institut 
dAstrophysique de Paris, which is gratefully acknowledged. 



REFERENCES 

Aguilar L.A., White S.D.M., 1985, ApJ. 295. 374 
Baugh CM., Cole S., Frenk C.S., [a.stro-ph/970311l| (submitted 
to ApJ) 

Binney J., Tremaine S., 1987, Galactic Dynamics. Princeton Uni- 
versity Press. 

Burstein D., Bender R., Faber S.M., Nolthenius R., 1997, AJ, 114, 
1365 

Chandrasekhar, S., 1943, ApJ, 97, 255 
Djorgovski S., Davis M., 1987, ApJ, 313, 59 

Dressier A., Lynden-Bell D., Burstein D., Davies R., Faber S.M., 

Terlevich R., Wegner G., 1987, ApJ, 313, 42 
Dressier A., Oemler A., Couch W.G., Smail I., EUis R.S., Barger 

A., B utcher H.. Poggiant i B.M., Sharpies R.M., 1997, ApJ, in 



press, istro-ph/9707232 



Geiger, B., & Schneider, P. 1997, MNRAS, submitted 

Ghigna S., Moore B., Governato P., Lake G., Quinn T., Stadel 

J., in preparation 
Gunn, J.E., Gott, J.R. Ill, 1972, ApJ, 176, 1 
Lavery R.J., Henry J. P., 1988, ApJ, 330, 596 

Kauffmann G., White S.D.M., Guiderdoni B., 1993, MNRAS, 264, 
201 

King I.R., 1962, AJ, 67, 471 



Klypin, A., Gottlober, S., Kravtsov, A. V. 1997, istro-ph/9708191 
Lacey C.G., Cole S., 1993, MNRAS, 262, 627 
Mamon, G.A., 1987, ApJ, 321, 622 

Mamon, G.A., 1993, in Combes P., Athanassoula E., eds., Af-body 

Problems and Gravitational Dynamics, 188. 
Moore B., Katz N., Lake G., 1996, ApJ, 457, 455 
Moore B., Katz N., Lake G., Dressier A., Oemler A.O.Jr., 1996, 

Nature, 379, 613 
Navarro J.F., Frenk C.S., White S.D.M., 1995, MNRAS, 275, 56 
Navarro J.F., Frenk C.S., White S.D.M., 1996, ApJ, 462, 563 
Oemler A., Dressier A., Butcher H.R., 1997, ApJ, 474, 561 
Spitzer L., 1958, ApJ, 127, 17 

Summers P. J., Davis M., Evrard A.E., 1995, ApJ, 454, 1 
Tormen G., 1997, MNRAS, 290, 411 

Tormen G., Bouchet F.R., White S.D.M., 1997, MNRAS, 286, 
865 

Weinberg M.D., AJ, 108, 1403 
White S.D.M., 1976, MNRAS, 174, 19 
White S.D.M., Rees M.J., 1978, MNRAS, 183, 341 
White D.A., Fabian, A.C., Forman, W., Jones, C, Stern, C, 1991, 
ApJ, 375, 35. 



© 0000 RAS, MNRAS 000, 000-000 



